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The measured transverse momentum spectra and HBT correlations of bulk (i.e. low pt) matter 
can be well explained by assuming that the soft sector of particles produced in ultrarelativistic 
heavy-ion collisions is (approximately) thermalized and undergoes collective accelerated expansion 
in both longitudinal and transverse direction. However, this implies that bulk matter will have a 
non-vanishing flow component transverse to the trajectory of a high pr partonic jets. In general, 
this will increase the energy loss experienced by the jet parton and modify the shape of the jet cone. 
In this paper, we present a systematic study of the magnitude of the additional energy loss induced 
by flow under realistic assumptions for the medium evolution. We argue that a perturbative QGP 
description may be sufficient for the measured Raa if flow during the medium evolution is taken 
into account properly. 

PACS numbers: 25.75.-q 



I. INTRODUCTION 

Energy loss of a high px 'hard' parton travelling through 
low pt 'soft' matter has long been recognized as a promis- 
ing tool to study the initial high-density phases of ultra- 
relativistic heavy-ion collisions (URHIC) [1 0, H 0, H H 
In pi-ITcH- it has been suggested that a flow compo- 
nent transverse to the high px parton trajectory would 
lead to increased energy loss as compared to the one in 
a static medium. 

Descriptions of measured low px spectra and HBT cor- 
relation radii for Au-Au collisions at RHIC energies in 
terms of a thermalized medium characterized by a col- 
lective expansion (flow) velocity field at thermal decou- 
pling are quite successful, see e.g. [llj. Experimental 
measurements of charge independent two particle corre- 
lations seem to indicate that the distribution of jet sec- 
ondaries is widened in longitudinal direction |l2j | . A pos- 
sible explanation could be that longitudinal flow leaves a 
characteristic imprint on the energy loss phenomenon. 
However, since jets propagate through the medium from 
the early onset of thermalization and lose energy mainly 
in the initial hot and dense phase, the knowledge of the 
flow profile at thermal decoupling is only of limited value. 



fnstead, one has to study a more generalized framework 
which allows to characterize the whole evolution of ther- 
malized matter. It is the aim of this paper to provide 
such a study in a framework which has been shown to 
lead to a good description of the experimental data of 
transverse mass spectra and HBT correlations 0] an d 
which is flexible enough to test also various other evolu- 
tion scenarios (which are at present not favoured by the 
data). In the present paper we focus on the energy loss 
of the leading jet particle and leave predictions for the 
distortion of the jet cone to a subsequent publication. 
II. THE FORMALISM 



The parton's energy loss depends on the position of its 
production point fo and the angular orientation of it's 
trajectory r. In order to determine the probability for a 
hard parton P(AE) to lose the energy AE while travers- 
ing the medium on it's trajectory, we make use of a scal- 
ing law 0] which allows to relate the dynamical sce- 
nario to a static equivalent one whose density has been 
rescaled. This probability distribution of the energy loss 
AE is given as a functional of the distribution u 4^ of 
gluons emitted into the jet cone |15|: 
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The explicit expressions for this quantity are different if 
one calculates under the assumption that multiple 
soft scattering processes or a few hard scattering pro- 
cesses are responsible for the energy loss. For both cases 
they can be found in 0] . Using these results we obtain 
P(AE) from the key quantity of jet energy loss, namely 
the local transport coefficient q(rj s ,r,T). The transport 



coefficient characterizes the squared average momentum 
transfer from the medium to the hard parton per unit 
pathlength. Since we consider a time-dependent inho- 
mogeneous medium, this quantity depends on spacetime 
rapidity rapidity rj s = \ ln(i + z)j(t — z), radius r and 
proper time r = \/t 2 — z 2 . For simplicity we focus on 
central collisions and assume azimuthal symmetry. The 
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transport coefficient is related to the energy density of 
the medium as q — ce 3 / 4 . In the case of an ideal quark- 
gluon plasma (QGP) one would expect c ss 2 [Ty] . It 
depends also on the jet trajectory r = fo + £n (produc- 
tion point ro and orientation ft) . In order to calculate the 
transport coefficient in the presence of flow, we follow the 
prescription suggested in [T(j and replace 
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where /3j_ is the spatial component of the flow field or- 
thogonal to the parton trajectory. T n±n± indicates the 
component of the energy-momentum tensor, where n±_ is 
orthogonal to the jet's trajectory. In the absence of any 
flow effects the original result q = ce 3 / 4 (p) is recovered. 
The transport coefficient enters the calculation of P(AE) 
via lo c and (qL) , which are the linearly line-averaged char- 
acteristic gluon energy: 
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and the time-averaged total transverse momentum 
squared: 
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respectively [l5j . 

The calculation of the transport coefficient and subse- 
quent u> c (ro, </>), (qL)(ro, <f>) and finally P(AE) requires 
a model of the fireball evolution which determines the lo- 
cal energy density and the flow profile. We discuss details 
of that model in the next section. In our calculation we 
average over all possible angles and production vertices, 
weighting the distribution of the jets with the nuclear 
overlap factor TaaO 3 ) = J dzp 2 (b,z) (for central colli- 
sions) with p the nuclear density as a function of impact 
parameter b and longitudinal coordinate z. 
This formalism sets the stage for the calculation of the 
nuclear modification factor. In the case of central colli- 
sions it is given by 



Raa{pt,v) = 



d 2 N AA /dp T dy 
T AA (0)d 2 a NN /dp T dy- 
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To obtain Raa we use the formalism and averaging pro- 
cedure as described above and calculate the inclusive 
charged hadrons production in LO pQCD. This amounts 
to folding the average energy loss probability into the fac- 
torized expression for hadron production (for brevity we 
give schematical expressions, an explicit representation 
can be found in 0,0): 
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Here, fi/A(x,Q 2 ) denotes the distribution of the parton 
i inside the nucleus as a function of the parton momen- 
tum fraction x and the hard scale Q 2 of the scattering. 
Likewise, D v ^ 7T {z 1 (i 2 ^) denotes the fragmentation func- 
tion of a parton / into a hadron with the hadron taking 
the fraction z of the parton momentum and a fragmen- 
tation scale /j, 2 which should be a typical hadronic scale 
0(1 GeV). The expressions for the hard pQCD cross 
sections aij^f+k can e.g. be found in I19| . W e use the 
CTEQ6 parton distribution functions [2£j, |21j for the pp 
reference, the NPDF set [2^ for production in nuclear 
collisions and the KKP fragmentation functions [23| . 
Unless stated otherwise, we assume for the strong cou- 
pling a s — 0.3 in the calculation of P(AE) in following. 



III. THE FIREBALL EVOLUTION MODEL 

The fireball evolution enters this framework in the shape 
of e(r] s ,r, r) and the flow profile u M (?y s , r, t). For the de- 
scription of the evolution, we base our investigation on 
the formalism outlined in |13|. 

The main assumption for the model is that an equili- 
brated system is formed a short time To after the on- 
set of the collision. Furthermore, we assume that this 
thermal fireball subsequently expands isentropically un- 
til the mean free path of particles exceeds (at a time r/) 
the dimensions of the system and particles move without 
significant interaction to the detector. 
For the entropy density at a given proper time we make 
the ansatz 



s(r,r) s ,r) = NR(r,r) ■ H(r) s ,r) 



(9) 



with t the proper time as measured in a frame co-moving 
with a given volume element and R(t,t), H{t] s ,t) two 
functions describing the shape of the distribution and N a 
normalization factor. We use Woods-Saxon distributions 



R(r, t) = 1/ [I + cxp 
H(t] s ,t) = 1/ (l + cxp 
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to describe the shapes for a given r. Thus, the ingredi- 
ents of the model are the skin thickness parameters d ws 
and 77 ws and the parametrizations of the expansion of the 
spatial extensions R c (t), H c {t) as a function of proper 
time. For a radially non-relativistic expansion and con- 
stant acceleration we find R c (t) = R + ^-t 2 . H c (t) is 
obtained by integrating forward in r a trajectory origi- 
nating from the collision center which is characterized by 
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a rapidity T) c (t) = rjo + a^r with r\ c = atanh v c z where 
v c z is the longitudinal velocity of that trajectory. Since 
the relation between proper time as measured in the co- 
moving frame and lab time is determined by the rapidity 
at a given time, the resulting integral is solved numeri- 
cally (see [24| for details). Rq is determined in overlap 
calculations using Glauber theory, the initial size of the 
rapidity interval occupied by the fireball matter. 770 is a 
free parameter and we choose to use the transverse ve- 
locity = a±Tf and rapidity at decoupling proper time 
V* = Vo + o-jjTf as parameters. Thus, specifying rio, rjf, Vj_ 
and Tf sets the essential scales of the spacetime evolution 
and (i W s and 7/ws specify the detailed distribution of en- 
tropy density. 

For transverse flow we assume a linear relation between 
radius r and transverse rapidity p — atanh uj_(r) = 
r/R c (r) ■ Pc(t) with p c (r) = atanh a±r. In the following 
investigation we also test a flow profile p ~ r 2 . 
We allow for the possibility of accelerated longitudinal 
expansion which in general implies r\ ^ r\ s [24j| . Here, 
71 = i In Pa+Pz denotes the longitudinal momentum ra- 
pidity of a given volume element. We can parametrize 
this mismatch between spacetime and momentum rapid- 
ity as a local Ar/ = r\ — r\ s which is a function of r and 

Vs- 

IV. RESULTS 

In the following we calculate in the multiple soft 
scattering approximation. The largest remaining uncer- 
tainty is the detailed choice of the parameter c in the rela- 
tion q = ce(T n±n± ) connecting transport coefficient and 
local energy density. As the highest temperature reached 
in the model evolution is only approximately two times 
the phase transition temperature Tq there is no a priori 
reason to assume that the ideal QGP value c « 2 would 
be realized. Hence, unless stated otherwise, we adjust c 
such that the scenario without flow reproduces the large 
Pt region (we do not include any intermediate px physics 
such as the Cronin enhancement in this study) and dis- 
cuss changes induced by flow relative to this baseline. 
As in 0| we find that Raa stays rather flat as a function 
of pt out to more than 50 GeV. If the flattening which 
may be seen in the data can be trusted, a value of c = 4—6 
would lead to agreement with the data by the STAR 
collaboration for 5% central 200 AGeV Au-Au collisions 
[25| (see Fig. QJ. Since taking into account the effect 
of flow transverse to the jet axis is expected to increase 
energy loss, we choose the value c = 4 as a baseline for 
the following investigation. 

Taking the transverse flow field into account we indeed 
observe increased supression, i.e. a reduction of Raa- 
Its magnitude shows some dependence on how the flow 
field depends on the radial position, we investigate p ~ 
r and p ~ r 2 . In both cases the maximum transverse 
velocity is adjusted such that the measured 7r~, K~ and 
p transverse mass spectra are reproduced. The emerging 
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FIG. 1: Calculated Raa as a function of transverse momen- 
tum for different values of the proportionality constant c be- 
tween transport coefficient q and energy density e 3,/4 com- 
pared with STAR data [25f . all without taking the effect of 
flow to the energy loss into account. 

trend is that shuffling more flow towards the fireball edge 
leads to a greater reduction of Raa- The reasons are 
rather involved: In order to experience transverse flow 
at all, the jet production vertex has to be away from 
the center of the transverse plane and the hard parton 
needs to propagate transverse to the local radial direction 
e r . On the other hand, hard partons produced close to 
the surface tend to escape before a significant amount 
of transverse flow could be generated. In the end, hard 
partons still in the medium at relatively late times will be 
close to the surface and experience additional energy loss 
if the flow profile leads to larger flow for the periphery. 
The main reason that the influence of transverse flow 
is much less pronounced than in the findings of [Tg| is 
that in a realistic evolution scenario the transverse flow 
field does not build up instantaneously but rather devel- 
ops gradually over time — as jet energy loss is a phe- 
nomenon predominantly sensitive to earlier times when 
the medium energy density is large, the spacetime re- 
gion of the evolution where transverse flow is small is 
probed with Raa- We demonstrate this in Fig. [5] where 
we vary the amount of pre-equilibrium transverse flow 
v T at the fireball edge. (We take care that the final 
transverse flow agrees with the measured m-r spectra by 
reducing the flow gained during the thermalized evolu- 
tion accordingly). The amount of primordial transverse 
flow is expected to be small, but our results demonstrate 
clearly that we recover the comparatively large sensitiv- 
ity seen in |10| when we assume that a strong flow field 
is present ab initio. 

However, just increasing the amount of primordial trans- 
verse expansion is hardly an appropiate comparison as 
in addition to the increased role of flow in the energy 
loss there is a rather trivial geometrical effect: The fast 
initial expansion of the fireball radius in the presence of 
primordial flow increases the average pathlength of hard 
partons in matter. In order to separate the two effects, 
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FIG. 2: Left panel: Calculated Raa without transverse flow (solid) and with a flow profile linear (dashed) and quadratic 
(dash-dotted) in radius as compared to STAR data pq] Right panel: Calculated Raa with linear flow profile for different values 
of pre-equilibrium flow velocity at the fireball edge Vj- . 



we show a calculation in which only the geometrical ef- 
fect was taken into account but where we neglected the 
flow effect on q. The results are shown in Fig. [3J 
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FIG. 3: Calculated Raa without primordial transverse flow 
Vt (solid black) and for 0.2 (blue) and 0.5 (red) showing 
just the geometrical effect of increased radius (dashed) and 
the full effect of geometry and increased energy loss (solid). 

In essence, for not too large primordial expansion veloc- 
ities, geometry accounts for about half of the observed 
effect. Note that a small primordial flow field cannot 
be ruled out by current experimental evidence: One of 
the observables most sensitive to the initial state is the 
emission of thermal photons, but using the framework 
outlined in 26] we observe that the increased radial flow 
tends to compensate the faster cooling by the increased 
expansion so that it is very difficult to distinguish the 
scenarios. 

In principle, in a non-boost invariant accelerated longi- 
tudinal expansion there can be a component of longitu- 
dinal flow transverse to the jet axis (this is not so in a 
Bjorken expansion where the hard parton is always co- 
moving with the surrounding matter in longitudinal di- 



rection). The magnitude of this component depends on 
the amount of longitudinal acceleration the system has 
undergone and the distance of the jet from midrapidity. 
In the present framework, it turns out that longitudinal 
flow leads to an additional 2% suppression at y = 0.5 
as compared to a situation in which only transverse flow 
is taken into account. This is small compared to other 
uncertainties, hence we do not discuss the effect of lon- 
gitudinal flow with regard to additional energy loss any 
further. 



V. THE RELATION BETWEEN ENERGY 
DENSITY AND TRANSPORT COEFFICIENT 

In |17| . an estimate within a Bjorken expansion model 
was made for the value of the parameter c relating the 
transport coefficient and energy density. The analysis 
found c > 8 ... 19 far from the perturbative value c ps 2 
valid for the ideal QGP. This is in fact consistently larger 
than our findings neglecting flow. The reason is that we 
do not use a Bjorken expansion. This leads to stronger 
initial compression of matter (and thus higher densities) 
as well as slower subsequent cooling in our model. If 
we assume a Bjorken expansion and neglect the effect of 
transverse flow on energy loss we find good agreement 
with the data for c = 10, quite consistent with the esti- 
mate in [l7|. 

However, if we take transverse flow into account in a 
moderately optimistic scenario, i.e. assuming a quadratic 
flow profile, a small value of primordial flow with v l T = 0.1 
and furthermore increase a s to 0.45 instead of 0.3 then we 
find that c = 2 in fact gives a fair description of the data. 
Thus, there may not be a huge quantitative discrepancy 
between pQCD predictions and the measured energy loss. 
This is shown in Fig. 0] 
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FIG. 4: Calculated Raa in a Bjorken expansion and without 
the effect of transverse flow on energy loss, requiring c = 10 
(solid red) and in the best fit evolution described in [13J, as- 
suming a quadratic flow profile, a small primordial flow ve- 
locity = 0.1 and a s — 0.45. In this case, only c = 2 is 
required for agreement with the data. 



buildup of the transverse flow field over time. Jets need 
to be produced away from the center of the transverse 
distribution in order to see a component transverse to 
the jet direction, but if they are produced too close to 
the edge they are likely to escape before flow could build 
up. As a result, we find great sensitivity to the presence 
of a primordial flow field. 

We have also studied the role of a longitudinal flow com- 
ponent orthogonal to the outgoing jet. We observe that 
away from midrapidity additional energy loss is induced 
if the jet is not longitudinally co-moving with the ther- 
malized matter, however around midrapidity this effect 
turns out to be small compared to the influence of trans- 
verse flow. 

If we include the additional energy loss induced by trans- 
verse flow in a favourable (but still realistic) scenario, we 
find that the ideal QGP relation q sw 2e 3 / 4 gives a fair 
description of the data. Thus, deviations from the ideal 
QGP energy loss may not be large even in the tempera- 
ture range < 2Tc probed at RHIC energies. 



VI. SUMMARY 

We have investigated the effect of flow transverse to the 
jet direction on the energy loss of the leading parton. The 
net effect arises from an intricate interplay of several dif- 
ferent effects, among them the quadratic pathlength de- 
pendence of energy loss which tends to emphasize the role 
of late times, the expansion and cooling matter which di- 
lutes the system and emphasizes early times and gradual 
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